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A sterile neutrino witli a mass of a few IceV can play the role of a warm dark 
matter (DM). This can be realized in seesaw models with 3 left- and 3 right-handed 
neutrinos. It is possible to identify the keV neutrino to be one of the right-handed 
neutrinos leaving the other two to be much more heavier, the uSM model. We 
show that with this realization of keV neutrino DM, the model has an approximate 
Friedberg-Lee symmetry providing a natural explanation for the lightness of the 
right-handed neutrino. We also find that in this model the mixing parameters couple 
light and heavy neutrinos are strongly correlated, and can be large enough to have 
testable effects at the LHC for the two heavy right-handed neutrinos to be in the 
hundred- GeV range. 

PACS numbers: 14.60.Pq, 13.15. +g 
Introduction 

Cosmological and astrophysical studies show that there are dark matter (DM) in our 
universe. The DM contributes about 20% of the energy in our universe. The identity of DM 
is still not known. Many models have been proposed. Neutrino has long been considered 
to be one of the possible candidates. The left-handed neutrinos, the active neutrinos, have 
standard model (SM) weak interaction, and were in thermal equilibrium in the early universe. 
Active neutrinos with masses in the range of a few tens of eV to a few GeV would over close 
the universe and are therefore ruled out as DM. Active light neutrinos of mass less than 
a few tens of eV has problem with structure formation. Data constrain the sum of the 
three light masses must be less than an eV or so. An active neutrino is unlikely to play a 
significant role for DM. However, a right-handed neutrino with an appropriate mass and 
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a mixing with active neutrinos can play the role of DM. The right-handed neutrino does 
not have SM interactions and may not be in thermal equilibrium in the early universe. But 
in general they mix with left-handed neutrinos and may be produced in the early universe 
through oscillations of different types of neutrinos. With appropriate mass and mixing, a 
right-handed neutrino of mass a few keV can contribute the correct relic density of our 
universe. This uji DM belongs to the warm DM category. 

Right-handed neutrinos can be introduced in different ways. If right handed neutrinos 
have large Majorana masses, they can play a very important role to explain why the light 
neutrinos have very small masses via the seesaw mechanism A priori, the right-handed 
neutrino mass scale is not known which can be as high as the Planck scale or as low as 1 



eV leading to many interesting consequences 



5|. We will study some implications of a 



seesaw model with a keV mass right-handed neutrino playing the role of DM companied by 



two heavy neutrinos. A model, called z/MSM to realize this has been proposed in 



3|. It is 



a seesaw model jl| in which there are 3 left- and 3 right- handed neutrinos, the 3-1-3 model, 
with one of the right-handed neutrinos (z/r, ) having a mass of a few keV. This is a minimal 

n n 

model of this type. In the model proposed in Refs. [2|, |3| the heavy neutrinos iyR2,3 have 
masses around 1 — 10 GeV. In the present work we take a different approach to have 1^/223 
mass scale to be in the hundred-GeV range and to study some implications for LHC physics. 
We refer this model as i/SM. 

In order for the keV right-handed neutrino to play the role of a warm DM, the parameters 
in this model are constrained. We show that the model has an approximate Friedberg- 
Lee symmetry providing an natural explanation for the lightness of the keV right-handed 
neutrino. We also find that in this model the mixing parameters which couple light and heavy 
neutrinos are strongly correlated, and can be large enough to have testable consequences at 
the LHC. 

Dark matter in z/SM 

Now we briefly describe the warm DM in z/SM. Let us indicate the keV scale right-handed 
neutrino which plays the role of the DM as z/r^, and the other two right-handed neutrinos 
which have masses in the multi-GeV or higher region as z/jijjg. The other particles relevant 
to our discussions are the three generations of left-handed lepton doublets = (z/^j, , 
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and the Higgs doublet H = [H'^, H )-^. The Lagrangian responsible to neutrino masses is 

C = -UrMv^r - LlYHur + H.C. , (1) 

where Ur is the charge conjugate of ur. M and Y are 3x3 matrices. M is the Majorana 
mass matrix of ur and is symmetric. For convenience we will work in the basis where M is 
diagonal, that is, M = diag{Mi, M2, M3}. 

In seesaw models one usually assumes that the right-handed neutrinos are super-heavy 
in the 10^^ ~ 10^^ GeV range, so that the masses of light left-handed neutrinos are strongly 



suppressed. But the scale need not to 
few hundreds of GeV to a TeV. In Ref. 



De so high. In fact the scale can be as low as a 
2| it was suggested that the lightest right-handed 
neutrino ur^ can even have a mass a few keV and play the role of a warm DM. 

After the electro-weak symmetry breaking, that is, the Higgs develops a non-zero vacuum 
expectation value {H) = {v, 0)'^, the neutrino mass matrix in the basis (i^l, lyR)'^ is given by 

(2) 

Diagonalization of the above seesaw mass matrix leads to the reduced mass matrix rriiy for 
light neutrinos, 




m,, = —V 



The above seesaw mass matrix also leads to a non-zero mixing matrix between the left- and 
right- handed neutrinos given by 

R = {Ru) = Yv{M*)-\ (4) 

where / runs over the light neutrinos, i^e.^i.r and i runs over, ur^ ^ 3 • 

For keV scale ur-^, the above mixing can cause oscillation, in the early universe, between 
the right- and left- handed neutrinos and lead to a finite energy density from the keV right- 



handed neutrino 
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Eq. is for production without resonance. For production with resonance, sufficient 
amount of keV scale uri warm DM can be produced for a much smaller mixing [81]. 
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It is interesting that uj^^ of a few keV can provide the right amount of warm DM density 
if it is stable enough. The keV scale Ufi^ can only decay into light active neutrinos at the 
tree-level through mixing of left- and right-handed neutrinos and exchanging Z and Higgs. 
At loop level it can also decay into photon through exchanging W and charged leptons. The 
tree-level decay modes of z/^^ are: i^r^ u + 2u and Ufi^ 2v + v. The lifetime of vr-^ is 
estimated to be {2! 

One can see that the lifetime of is much larger than the age of the universe ~ lO^^s for 
mass Mri of order a keV and \Rii\^ ~ 10~^. A keV scale vpii has a long lifetime which 
is allowed for a DM candidate. 

Constraints on the mass and the mixing of vri as warm DM candidate come from struc- 
ture formation {2!, phase space density jsl, X-ray lines of i^a{i^a) + 7 decay 10-12| 
and Lyman-a forest data 13||. Astrophysical implications of keV scale ur, such as effect on 
supernova explosion and re-ionization etc., have been explored in Refs. IJ, ll5| . A review 
on related subjects can be found in Ref. [161]. Considerations on structure formation, phase 
space density and Lyman-a forest data give lower bound on the mass of z/^j;^ warm DM. 
Observations of X-ray lines of z/^^ decay give upper bound on the mass and mixing of z^r^ . 
These constraints can be very strong if t/p, warm DM accounts for all the DM density in the 



universe. However, a recent analysis 



131 ] shows that these constraints can be significantly 



released for < 0.4 Qdm- For example, z/r^ with mass in the range 1 keV < Mi < 5 keV 
and \Rii\'^ ~ 10"^ — 10~® is allowed by all these constraints if ^uri = 0.2 Qdm- We will 
use 1 keV < Mi < 5 keV and \Rii\'^ ^ 10~® for discussions later. 

The other two right-handed neutrinos, uji^ 3 , can have large masses and generate the 
light active neutrino masses and mixing through the seesaw mechanism. This model can 
correctly account for the masses and mixing of light neutrinos by making appropriate choice 
of Yukawa couplings 2|] . 

Implication of the keV neutrino for the Yukawa couplings 



In general the Yukawa coupling matrix Y can be parameterized in the following form 17 \ 



Y = -U{miy/^0{M*y/\ (7) 

V 
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where = diag{mj^^e*'^^/^, m2^^e*'^^''^, mg^^e*'^^''^} where real numbers mj(z = 1,2,3) are 
the masses of three hght neutrinos and </)j(z = 1, 2, 3) are three Majorana phases. (M*)^/^ = 
diag{(M*)i/^ (M*)i/2^ (^3*)^^^}- U is the neutrino mixing matrix observed in experiments, 
•fhi, is the diagonahzed neutrino mass matrix, O is a complex orthogonal matrix: 0^0 = 
00^ = 1. 

Using Eq. ([7]), the mixing matrix R can be written as 



R = U{mlf/^0{M 



r*\-l/2 



(8) 



The orthogonal matrix O can have determinant det{0) = ±1. Orthogonal matrix O with 
det(O) = +1 can be parameterized as the products of three matrices: 



O = O23O13O12, 



(9) 



where Oij is a 3 x 3 complex orthogonal matrix with determinant +1. It can be expressed 
using complex numbers 9ij, e.g. 



^ cos6'i2 sin ^12 0^ 



O 



12 



V 



(10) 



sin 6^12 cos 6^12 

V 

Orthogonal matrix O with det(O) = — 1 differs from an orthogonal matrix expressed in 
Eq. ([9]) by a transformation using matrices such as diag{l, 1,-1}, diag{l, —1,1}, etc. These 
factor —1 can be absorbed into the phase factors in (m*)^/^. That is, we can choose 

< 0i < 47r (11) 

instead of in the range [0, 27r] to account for the possibility that det(O) can be —1. 

To further understand implications of i^SM we rewrite the seesaw mass formula, Eq. (IHl), 
as the following 



1,2,3 



(12) 



where 



(13) 
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Eq. (jl]) has been used in obtaining Eq. (fT5]) . Si is the contribution of to the hght neutrino 
masses. Using the constraint l-R/iP ~ 10"*^, we find for 1 keV ^ Mi < 5 keV 2! 



We introduce 



Y' 
R' 



Using Eq. (fTS!) we get 



The condition Eq. ([T 



a6 



^10"' eV. 

= U^R = mll^O{M*)-^l'', 
^ U^S^U = -v'Y;:Y;*Mr\ 



(14) 



is re-expressed as 

\iS[)a,\<10-' eV 



The experimental constraints on the neutrino masses, at 2a level 

7.25xlO~^eV^ < Am^^ < S.llxlQ-^eV^ 
2.18xlO~^eV^ < \Aml^\ < 2.64 x 10"^eV^ 



are 



(15) 

(16) 

(17) 

(18) 
(19) 



^^mj < 1.2 eV. For normal hierarchy(NH) of neutrino masses, m2 ~ a/ Am|i ^ 0.9 x 10 ^ 



eV and ma ^ \/TA 



m 



311 



That is 



0.05 eV, Eq. ([HD says \{S[] 



|Oai|<l, a = 2, 3 



"^al^^all ^ for O = 2, 3. 



(20) 



From this condition we find that the matrix O is approximated as 



O 











COS 9 sin 6 
— sin 6* cos 9 



(21) 



For inverted hierarchy(IH) of neutrino masses, mi ^ a/| Am|i| 0.05 eV and m2 
Am2^ ^ 0.05 eV, ^ says |(5;)aa| = ^a|0*f | < m^ for a = 1, 2. That is 



m 
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|Oai|<l, a = 1,2 



(22) 
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From this condition we find that the matrix O is approximated as 

^0 cos^ sin 6'^ 

-sin^ cos^ . (23) 
yl ) 

The above analysis show that z/^j^ gives neghgible contribution to hght neutrino masses. 
This condition imphes that hght neutrino masses are dominated by contributions of two 
heavier right-handed neutrinos i'r2.3- The lightest left-handed neutrino has mass of order 
~ 10"^ eV. uSM can not reproduce degenerate mass pattern of light left-handed neutrinos. 
We also have a good idea about some of the properties of the Yukawa coupling matrix 
Y. We will show in the following that the resulting mass matrix exhibits an approximate 



Friedberg-Lee(FL) symmetry 19|, |20| 



The approximate Priedberg-Lee symmetry in uSNL 

A theory is said to have a FL symmetry when the Lagrangian of this theory is invariant 
under the transformation on a fermion field of the form q ^ q + e, where e is a space-time 
independent element of the Grassmann algebra, anti-commuting with the fermionic field 
operators q. We explain why the neutrino mass matrix in the uSM has an approximate FL 
symmetry in detail in the following. 

The general Yukawa matrix can be written as 



Y 



Y^i Yfj_2 Yfj_3 



(24) 



Using Eqs. ( 12T]) and (!23|) we find that for NH case, for Y part of the Y one has 

/ n n \ 



vY = U{ml) 




cos 6 sin 6 
— sin 9 cos 6 



a/2 



(25) 



where M = diag{M2, M3}. 
For IH case Y, one has 



^ cos^ sin^^ 



vY = U^mlf^ 



sin 9 cos 9 




(M*) 



*\l/2 



(26) 
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while the Y part of Y is given by 



/ (M*)l/2^ 



vY = U{ml) 



1/2 



v 








or vY = U{mlf/'^ 



J 



( \ 





(27) 



for NH or IH, respectively. 

Using result from Eq. fll7l) we can also write the light neutrino mass matrix as 



m,, —V 



'^Y*M^^Y'^ . 



(28) 



2l|. 



The above is the mass formula in the minimal seesaw model ^ 

^ 1 /2 1 /2 

Note that since |M2,3| ^ |Mi| and the Yw is suppressed by jM^^^ in comparison with 
elements of Y . In Eqs. (12T]) and ( 123|1 there are no large factors in the first columns to enhance 
the Yukawa couplings, Y is the leading term in Yukawa coupling, we have approximately: 



Y 



^0 y;2 




yo Yri 

The mass matrix M can be approximated as 



(29) 



M^diag{0,M2,M3}. 



(30) 



It is easy to check that with the above Y and M, the Lagrangian in Eq.([T]) is invariant 
under a PL transformation 



(31) 



If the PL symmetry is a global one, that is, e is space-time independent, it is easy to check 
that the kinetic term Ck = '^_r7/j(^<9^i^r) is also invariant under the transformation defined 
in Eq.(|3T]). 



In fact it has been shown 
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23l ] that imposing a global PL symmetry on a particular 



direction of right-handed neutrinos in the 3+3, it results in a massless right-handed neutrino, 
and the mass matrices is equivalent to the from of a 3+2 (two right-handed neutrinos) 
minimal seesaw model which predicts a light neutrino with a zero mass. In our case, uji-^ 
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is the corresponding massless right-handed neutrino, and the hghtest hght neutrino mass of 
order 10~^ eV corresponds to the zero mass one in the exact FL hmit. 

We conclude that the z/SM has an approximate FL symmetry. Small violation of this 
symmetry provides a natural explanation why one of the right-handed neutrino has much 
smaller mass than other two heavy right-handed neutrinos. We would like to comment that 
the 3-1-3 seesaw model is the minimal model which is consistent with light neutrino masses 
and can have a warm DM candidate in the framework of seesaw mechanism. We note that 
the FL symmetry we found is a consequence of the experimental constraints. Some other 



symmetries can be imposed to obtain z/SM 



24|. 



The possibility of a large mixing between light and heavy sectors in z/SM 

The model considered in Ref. (2I has very small Yukawa couplings. They assume the 
Yukawa couplings are of order ~ mJ'^M^J'^ jv. If this is always the case, then even the heavy 
right-handed neutrinos are light enough, a few hundred GeV, to be produced at the LHC, 
the small mixing between light and heavy neutrinos makes it impossible to be detected. 
Fortunately this is not necessarily true. There are other possibilities!^. 

It is known that the matrix elements of the complex orthogonal matrix O can be large if 
6' is a complex number: 



9 = X + iy 



(32) 



where x and y are two real numbers. In this case | cos6'| and | sin6'| can be enhanced by a 
large factor e'^' if \y\ is large. Elements in the Yukawa coupling Y matrix and the mixing 
R matrix can be enhanced by large elements of matrix O. Tiny light neutrino masses are 
reproduced with large Yukawa couplings in the seesaw formula through fine tuning. In the 
following we explain in more detail how a large mixing between light and heavy neutrinos 
can be obtained in z/SM. 

Consider the NH case. Eq. (12T]) can be written as 
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For \y\ ^ 1, the first term in the right-handed side of Eq. 
factor e'^l and is the leading term. 



(33) 



is enhanced by the large 
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Using Eqs. ([7]) and f l33p we find for |?/| ^ 1 the leading term in Yukawa coupling is 



1/2 



vY, 
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(34) 
(35) 



|yii| are small as explained in Eq. ( 127|1 . 
For IH case we find 
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(36) 



Again for \y\ 3> 1, the first term in the right-handed side of Eq. ( l36l) is the leading term. 
The leading terms in the Yukawa coupling are found to be 



vYi3 = ie^"-+l^l(±zf/ami/'e^'^i/2 + Ui2ml^\-"f''/^){M;y/\ 
From Eqs. dM]), (ES]), ([37]) and dM]) we find that 



(37) 
(38) 



Using Eq. 



Yi3{m;)-'/^ ^ ±zy,2(M*)-^/l 

Eq. fl39|) is re-expressed as 



(39) 



(40) 



Using Eqs. fl39|) and fHOjl . we see in Eq. f|T3|) that there is a strong cancelation between 5*2 
and 5*3, contributions of vr2 and Vfci to light neutrino masses. 

It has been shown in Ref. {26 1 that neutrinoless double decay gives strong constraint on 
the mixing and masse of sterile neutrino. Typically constraint for a single sterile neutrino 
is found to be |-ResP ^ 10"^ for Ms ~ 100 GeV. We note that this constraint does not 
apply to i'R2,3 of degenerate or quasi-degenerate masses. GeV scale i'r2,3 contribute to the 
neutrinoless double beta decay with the amplitude 



A = F{R^ 



■e2 



1 

M2 



R 



e3 



(41) 
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where F is a factor containing all other effects and M2,3 have been chosen real for convenience 
of later discussion. Eq. fHTj) can be rewritten as 



A = -^iRl2M2 + RlM,) + FRl,M,{^^ - ^). (42) 



According to Eq. ( I40l) and Eq. (1131) . the first term in Eq. (H2l) has a strong cancelation 
and is of order (F/M|)(m;^)ee and is negligible. The second term in Eq. (142|) is essential 
for constraining the z/SM. It's clear that if vr2,z are degenerate or quasi-degenerate, their 
contribution to neutrinoless double decay is greatly reduced. For example, if (AMfg/Ml) ~ 
10~^ where AMI2 = M| — M|, the rate of double beta decay is reduced by 10 orders 
of magnitude and the constraint becomes |-Re2.3p ^ 1- We note that a quasi-degeneracy 
of heavy neutrinos is required for producing baryogenesis using oscillation of right-handed 



28 



with resonance at the electro-weak scale. 



neutrinos js], Q or using thermal leptogenesis 

Matrix Y can be transformed using bi-unitary transformation to a diagonalized form: Y = 
diagjyi, 1/2, Z/s}- Large matrix elements in Eqs. ( !33l) and ( l36l) affect 7/2,3 which correspond to 
two linear combinations of quasi-degenerate i^r2,3- Note that the matrix O does not change 
the determinant of Y. Hence one of 7/2,3, which we take to be 7/3, is enhanced and another 
one, 7/2; is suppressed. The fact that only one of 7/1,2,3 is enhanced can be seen clearly in 
Eqs. (l33|) and (136|) . that is the leading part of the O matrix, enhanced by e'^', is of rank 
one. Hence one can make 7/2 ~ 10^^ and 7/3 large. This means one of the linear combination 
of quasi-degenerate i'r2,3 can never be in thermal equilibrium before the freeze-out of the 
sphaleron transition which is necessary for producing baryogenesis. Whether with large 
elements in O realistic baryongenesis can be obtained is an interesting question to study. 
Detailed study on this subject is out of the scope of the present article and will be studied 
elsewhere. Here we are interested to see if there are parameter spaces in which experiments 
at the LHC can probe. 

We summarize some of the interesting properties in the following: 

• The Yukawa couplings Yi2,i3 are enhanced by large e'^' if \y\ is large. 

• Only one of the Majorana phases (pi is observable. In both NH and IH cases we can 
set 01 = 03 = and keep 02. For convenience we write 

(m:)V2 = diag{mi/^mf e'*,mf }, (43) 
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where $ = — </)2/2. $ can be chosen in the range [0, 27r] as discussed for Eq. fITT]) . 

• None of the Majorana phases in Mi is observable since one can always rotate away 
phases of Mj. 

• A CP violating phase =i=x in the factor e^*^ is observable. 

• For GeV scale degenerate 2/7^2,3 the constraint from the neutrinoless double beta decay 
on the mixing R^i is greatly reduced. 

The right-handed heavy neutrinos do not have SM gauge interaction. If they do not mix 
with left-handed neutrinos or the mixing is extremely small, it is not possible to produce 
them and study their properties. As we have shown in the above that it is possible to 
have large mixing. The production of heavy neutrinos at the LHC becomes possible. Using 
constraints from electro- weak precision test, it has been shown model-independently that 
the heavy neutrinos can be produced and studied at the LHC up to 400 GeV for 2cr with 
100 fb~^ of integrated luminosity and mixing \Rii\^ < 10"'^ {25 1. The mixing of vr-^ with 
left-handed neutrinos are too small to have direct laboratory observable effects. But vr^^ 
may be produced and their properties may be studied. We note that \Rei\^ ^ 10^^ can be 
reached in z/SM, consistent with constraint from double beta decay experiments, for quasi- 
degenerate vr2,z which have AMI2/MI ^ 10~^. In the later discussion of the decay pattern 
of z^_R,2,3 we will concentrate on the degenerate or quasi-degenerate case. 

Correlation of vr^ 3 and light neutrino mass hierarchy in z/SM 

As seen from Eqs. ([7]) and ([H]), for a given form of O, one can establish the connection 
between the Yukawa coupling F, the mixing matrix i?, and the properties of the light 
neutrinos {U and masses). The leptonic mixing matrix U is usually written as 



U 



^ C12C13 C13S12 e '^si3 ^ 

~Ci2Si3S23e'^ — C23S12 C12C23 — e"'si2Si3S23 C13S23 
— e^''ci2C23Si3 — £235125136"^ — C12S23 C13C23 



(44) 



y ■'12'^23 ~ e C12C23613 —(-23 ■'12 ■^136 — l-12^'23 '-13'-23 J 

where Sjj = sin6'jj, Cij = cos 6ij, < 6ij < 7r/2 and < 6 < 27r. The phase 6 is the Dirac 
CP phase. The experimental constraints on the neutrino masses are listed in Eq. fITI?]) and 
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FIG. 1: |i?i2p-^^2/100 GeV for NH (left) and IH (right), assuming vanishing Majorana phase and 
< X < 27r, -20 < y < 20. 



the constraints on the mixing parameters, at 2a level 18l|, are 



0.27 < sin^^is < 0.35, 

0.39 < sin2^23 < 0.63, (45) 
sin2^i3 < 0.040. 

Using the above experimental constraints, one can explore the allowed values for the 
Riiil = e, ^,T;i = 2, 3) couplings and the heavy masses. In Fig. [1] we plot the allowed values 
for the normalized couplings |-R;2p^2/100 GeV of each lepton flavor for each spectrum, the 
NH (left panel) and the IH (right panel), assuming vanishing Majorana phase. We see that 
the couplings of different lepton flavors have strong correlation as expected. 

Moreover, as mixing elements of matrix R govern processes involving right-handed neu- 
trinos, most stringent model-independent data from precision electro-weak measurements 
and low-energy lepton-number violating processes can constrain unknown parameters in the 
complex orthogonal matrix O and the Majorana phase We show |-Ri2p as a function 
of the imaginary part y in parameter 6* of O in Fig. O One can see that |-Ri2p can reach 
~ 10~^, the upper bound from precision measurement. Constraint from neutrinoless double 
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FIG. 2: |i?;2p versus parameter y in matrix O for NH (left) and IH (right), assuming vanishing 
Majorana phase and < x < 27r, M2 = 100 GeV. The horizontal lines represent the current exper- 
imental bounds from precision electro-weak measurements 



29( 1 and searches of neutrinoless double 
beta-decay, the dashed and dotted lines, assuming AM32/M2 to be 10~^ and 10~^, respectively . 

beta decay is weaker for quasi-degenerate i^_r,2,3- As noted previously, heavy neutrinos with 
l-R/jp ~ 10"^ can be produced and studied at the LHC 

The partial decay widths of the heavy Majorana neutrinos are given by {25 1 



9 in|2ji#3/'i ,. \2 



\RuVMt{l-^^. 



iW 



9 I 7-) |2 ji^ /I ,. \2 



327r 
9' 



iW 



9' 
327rc^ 



(46) 

(47) 
(48) 

(49) 



where /ijj = Mj/Mf. If ujn is heavier than the Higgs boson one has the additional 
channels 

,2 
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Ru\'Mf{l~fiiH) 



(50) 



As discussed above, the lepton-fiavor contents of decays will be different in each neutrino 
spectrum. In order to search for the events with best reconstruction, we will only consider 
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FIG. 3: Branching fractions of degenerate neutrinos Ufi2 — > l'^W~ {l~W^) (/ = e,fi,T) for NH 
(left) and IH (right) versus heavy neutrino mass with < x < 27r, y = 10 and Mh = 120 GeV, 
assuming vanishing Majorana phase. 

the decay to charged leptons plus a W^. In Fig. [3] we show the impact of the neutrino 
masses and mixing angles and parameters in O on the branching fractions of the heavy 
neutrino z//j2 decaying into e, fi, t lepton plus W boson, respectively, with the left panels 
for the NH and the right panels of the IH, assuming vanishing Majorana phase. We also 
plot the dependence of branching fractions on parameter y in Fig. |H For large values of 
heavy neutrino mass or y the branching fractions can differ by one order of magnitude in NH 
case BR{fi^W^), BR{t^W^) > BR{e^W^) and about a factor of few in the IH spectrum 
BR{e^W^) > BR{ij,^W^),BR{t^W^). 

In general, the decay rates depend on only one Majorana phase $ when the lightest 
neutrino mass vanishes in the NH or IH case. In Fig. |5l we show the dependence of z/r2 decay 
branching fractions on Majorana phase $ in NH and IH for y = 10. In NH the dominant 
channels swap from t^W^ when $ ^ 7r/2 to fi^W^ when $ ^ 37r/2 by a few times. In IH 
the dominant channels swap from e^W^ when $ ^ 7i/2 to ^^W^,t^W^ when $ ^ 37r/2 
by more than one order time of magnitude. Moreover, it is important to note that the curves 
of branching fractions corresponding to Majorana phase translate parallelly by a phase vr for 
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FIG. 4: Branching fractions of degenerate neutrinos ujn l~^W~ {l~W~^) (/ = e, fi, r) for NH (left) 
and IH (right) versus parameter y with < x <2tt, M2 = 300 GeV and Mh = 120 GeV, assuming 
vanishing Majorana phase. 

—y case. This qualitative change can be made use of extracting the value of the Majorana 
phase $ and parameter y. 

If nature indeed uses low scale heavy neutrinos, of order one hundred GeV, with large 
mixing to light neutrinos, they may be produced at the LHC The z/SM can be tested by 
studying correlations between the decay patterns of heavy neutrinos and light neutrino mass 
hierarchies. 

Conclusions 

In summary we have studied the symmetry and some phenomenologies of the z/SM. 
The minimal model of this type is a 3 + 3 seesaw model with one of the right-handed 
neutrino to be light to play the role of the keV warm DM. We found that the mass and 
mixing parameters of the DM right-handed neutrino, constrained by experimental data, lie 
in the range that there is an approximate FL symmetry in the Lagrangian of the z/SM. 
The masses of the light active neutrino mass hierarchy, predicted in i/SM, can be IH or 
NH, but can not be quasi-degenerate. The seesaw masses of the light active neutrinos 
are dominated by contributions of two multi-GeV right-handed neutrinos vr^^^ and can be 
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FIG. 5: Branching fractions of z//j2 l^W {I W~^) (/ = e, /u,r) versus Majorana phase ^ for NH 
(left) and IH (right) when M2 = 300 GeV, Mh = 120 GeV and < 2; < 27r, y = 10. 

approximated by 3 + 2 seesaw formula. We should emphasize that the approximate FL 
symmetry discussed in the present work is not only valid for the range of i'r2,3 masses 
discussed, it is also valid for seesaw model with much larger or smaller masses of u^i^ 3 , 
as long as the keV scale dark matter constrained by observational data is included in the 
seesaw model. The keV scale ujn has a number of astrophysical implications, such as effect on 
supernova explosion |l^ and on the re- ionization 15 1. The phenomenology of the two multi- 



GeV scale right-handed neutrinos is similar to the phenomenology of 3 + 2 seesaw model. 
The lightest light neutrino has a mass of order 10~^ eV. We found that in z/SM the Yukawa 
couplings can be large. In particular we found that double beta decay experiment does 
not give strong constraint on the mixings and masses of two multi-GeV scale right-handed 
neutrinos if they are degenerate or quasi-degenerate. The Yukawa couplings can be large 
enough to be tested in LHC experiments. We also found that there are strong correlations 
between the couplings of the two heavy neutrinos 3 light neutrino mass hierarchy. 
The decay patterns of these heavy neutrinos sensitively depend on the Majorana phase. The 
decay patterns of the right-handed neutrinos u^^ 3 can be used to extract information of the 
mass hierarchy and Majorana phase of the light active neutrinos. 
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Note added: After the submission of the present article, recent analysis on X-ray obser- 
vations of local dwarf Willian 1 show evidence that the z/r^ dark matter may have mass 
around 5 keV with mixing |-RnP ~ 10~^ 3o|. Another analysis on the X-ray observation 
of the galactic center suggests that v^i dark matter has mass around 17 keV with mixing 



~ 10~^^ 3l|. It's easy to see that Eq. f|T^ is valid for these two groups of parame- 
ters of vr^ and the approximate Friedberg-Lee symmetry discussed in this article applies to 
models with these ranges of parameters space. 
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